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In general, different countries and communities respond to epidemics in accordance with their own control 
plans and protocols. However, owing to global human migration and mobility, strategic planning for 
epidemic control measures through the collaboration of relevant public health administrations is gaining 
importance for mitigating and containing large-scale epidemics. Here, we present a framework to evaluate 
the effectiveness of random (non-strategic) and targeted (strategic) epidemic interventions for spatially 
separated patches in metapopulation models. For a random intervention, we analytically derive the critical 
fraction of patches that receive epidemic interventions, above which epidemics are successfully contained. 
The analysis shows that the heterogeneity of patch connectivity makes it difficult to contain epidemics under 
the random intervention. We demonstrate that, particularly in such heterogeneously connected networks, 
targeted interventions are considerably effective compared to the random intervention. Our framework is 
useful for identifying the target areas where epidemic control measures should be focused. 

Modern transportation and traffic networks have enabled people to travel more frequently and globally 
than before, and therefore, epidemic outbreaks in local regions can potentially threaten international 
public health. The diversity and complexity of human mobility patterns complicate the planning of 
effective control measures against global epidemics at a local level. Therefore, epidemic control planning through 
international and inter-community cooperation is increasingly gaining importance to mitigate the threat of 
public health emergency caused by global epidemics 13 . The challenges to improve the detection and containment 
of epidemics with international importance have been addressed in international surveillance and response 
programs 4 , such as Global Outbreak Alert and Response Network (GOARN) 5 supported by the World Health 
Organization (WHO). The collaboration networks of public health administrations and research institutes are 
expected to contribute to the planning of strategies and policies for better epidemic control measures. An 
assessment of the impact of potential control strategies is a critical issue for policymakers to implement them. 

A metapopulation model is a useful mathematical framework to understand the dynamics of epidemic spread- 
ing in spatially separated patches, corresponding to countries, cities, or local communities, which are connected 
by pathways for human migration and population mobility 6 9 . Data analyses revealed the heterogeneous struc- 
tures of transportation networks such as airline networks 10-12 , railway networks 13,14 , and commuting networks 15,16 , 
and therefore, metapopulation models with complex patch connectivity have been intensively studied 15 29 . A 
major concern in these studies has been the global epidemic threshold, which determines if a global epidemic 
occurs or not when a small number of infected individuals invade a local patch 16,19,20,23-25 . The metapopulation 
framework has also been widely employed to investigate the efficacy of containment measures such as travel 
restrictions 19,21,28,29 , vaccinations 17 , and therapeutic use of antiviral drugs 21 . Colizza et al. 21 demonstrated that the 
redistribution of the stockpiles of antiviral drugs between prepared and unprepared countries is effective under 
limited antiviral supplies. However, such a cooperative strategy for controlling the infection of subpopulations in 
local patches has yet to be fully examined. 

Here, we present a metapopulation framework to evaluate the effectiveness of strategic interventions for 
epidemic control as compared to a non-strategic intervention. We consider a public health intervention that 
significantly reduces the infection risk of a subpopulation in a local patch through various countermeasures, e.g., 
the surveillance of cases and health indicators, delivery of vaccines, distribution of antiviral drugs, and promotion 
of hand hygiene. Public health interventions are provided to a fraction p of the patches, which we refer to as the 
low-risk patches* 0 . We assume that the infection rate in the low-risk patches is lower than that in the other patches 
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without interventions. An increase in the ratio p of the low-risk 
patches generally leads to a decrease in the prevalence of infection. 
Under certain conditions, there exists a critical ratio p c of the low-risk 
patches, above which the containment of epidemics is successful. 
Because the prevalence can differ depending on how to choose the 
low-risk patches, the critical ratio p c can also depend on the order in 
which low-risk patches are chosen as p is increased. If we appropri- 
ately select the low- risk patches, then it would be possible to decrease 
the critical ratio p c , implying a reduction in the total cost of public 
health interventions targeted at the local patches. Thus, we use the 
critical ratio p c as a measure for evaluating the effectiveness of public 
health interventions. First, we analytically derive the critical ratio for 
a random (non-strategic) intervention where low-risk patches are 
randomly chosen from the whole patches. Second, we numerically 
examine the effect of strategic interventions where the low-risk 
patches are chosen based on available information. The results show 
that strategic interventions are highly effective compared to the ran- 
dom intervention in metapopulation models with heterogeneous 
connectivity. Finally, numerical simulations performed using a 
Japanese airline network demonstrate that our framework is useful 
to specify the best intervention strategy and its target areas. 

Random and targeted interventions for controlling epidemics 
have been widely argued based on both contact networks of indivi- 
duals and the metapopulation models with networked patches. 
However, the targeted immunization in the contact networks corre- 
sponds to removal of the static nodes (i.e. the individuals) which do 
not have dynamics 31 35 . Therefore, the results on the immunization 
strategies based on the contact networks cannot straightforwardly 
apply to the intervention strategies at the metapopulation level. The 
targeted intervention in our metapopulation framework changes the 
transmission rate in the local patches. Such a model framework is not 
found in most of the relevant studies on metapopulation models 
where the transmission rates (or the local basic reproductive num- 
ber) are assumed to be identical. As far as we know, there is only one 
study taking into account patch differences in the local transmission 
rates, where the epidemic threshold is analyzed with respect to the 
mobility parameter 30 . The threshold with respect to the fraction of 
low-risk patches, which we focus on in the present study, gives a new 
perspective. This viewpoint enables to examine the strategies for 
epidemic control based on local interventions, in a framework dif- 
ferent from the previous ones. Although we employ a deterministic 
model without stochastic effects and assume a simple mobility pat- 
tern and simple infection dynamics in order to highlight the useful- 
ness of our framework and analysis, the method is applicable to more 
complicated metapopulation models incorporating stochastic fluc- 
tuations, realistic human mobility patterns 16,27 , and realistic human 
behavior 28 . 

Results 

Metapopulation models with epidemic intervention. We consider 
epidemic spreading in a population of individuals, separated into 
subpopulations in spatially discrete patches connected by migra- 
tion pathways. The model is represented by a metapopulation 
model consisting of N patches, in each of which susceptible and 
infected individuals are homogeneously mixed 19,24,25 . The epidemic 
dynamics in each local patch is described using a susceptible- 
infectious-susceptible (SIS) compartment model 36,37 : In patch i for 
i = 1, N, a susceptible individual interacting with an infected 
individual turns into an infected one at a transmission rate /J,- and an 
infected individual recovers to be susceptible at a recovery rate fa. 
Note that, although we use the SIS model for simplicity in this paper, 
our framework is also applicable to other compartment models of 
epidemic spreading. Here we assume non-limited transmission for 
simplicity, in which the total number of contacts of the individuals 
within a patch increases with an increase in the population density of 
the patch 19,24,25 (see Supplementary Information for the theoretical 



analysis of the case with limited transmission, in which the contact 
rate is constant independently of the population density). If patch i is 
isolated, the epidemic threshold is given by the local basic 

reproduction number defined as Rq = ft J fij for the SIS model. In 
metapopulation models, however, the global invasion threshold 
depends on the other factors such as population density, patch 
connectivity, and mobility rates 19,24,25 . The total density of the 
individuals is represented by p. The connectivity of the patches is 
characterized by the degree fc, for i = 1, . . ., N, which represents the 
number of migration pathways for patch i. The mobility rates of 
susceptible and infected individuals are represented by D s and D h 
respectively. Each individual in a patch randomly chooses one of the 
pathways departing from the patch and migrates to another patch 
with the mobility rate. The dynamics of epidemic spreading in the 
metapopulation model is divided into two processes: transmission or 
recovery within local patches (the reaction process) and migration of 
individuals between patches (the diffusion process) 19,24,25 . If the 
reaction and diffusion processes occur simultaneously, then a 
suitable continuous-time formulation of the metapopulation 
dynamics can be derived 24,25 . We employ this continuous-time SIS 
metapopulation model. 

We incorporate public health interventions that lower the infec- 
tion risk in a local patch into the metapopulation model. The whole 
patches are divided into low-risk patches with interventions and 
high-risk patches without interventions. For simplicity, we assume 
that the recovery rate is common for all patches, i.e., ^, = fi, and the 
transmission rates are different between the two types of patches, i.e., 
P t = /? L < j.i for the low-risk patches and /?,■ = [3 n > [i for the high- 
risk patches. The local basic reproduction number for the low-risk 

patch is given by = fi L / ji< 1 and that for the high-risk patch by 

fH) 

_Rq — j8 H //i>l. Therefore, a disease-free state without infected 
individuals is stable in an isolated low-risk patch, whereas an 
endemic steady state with a certain number of infected individuals 
is stable in an isolated high-risk patch. In the metapopulation model 
where the low-risk and high-risk patches coexist, the steady state 
depends on the ratio between the two types of patches and their 
configurations. Figure 1 schematically illustrates the metapopulation 
model that consists of low-risk and high-risk patches. The ratio of the 
low-risk patches to the whole patches is represented byp with 0 s p 
£ 1. A large value of p indicates that the public health interventions 
are more costly in total. Let us suppose that the steady state is an 
endemic state forp = 0 and a disease-free state forp = 1. In such a 
case, it is expected that there exists a critical value p = p c , where a 
transition between the two states occurs. In this study, the critical 
low- risk patch ratio p c is used as a measure for evaluating the effec- 
tiveness of the interventions for epidemic control. 

Critical low-risk patch ratio for random intervention. We inves- 
tigate a random intervention where low-risk patches are randomly 
chosen from the whole patches. For the SIS metatpopulation model, 
the epidemic size is evaluated by the prevalence of infection Pi/p, 
where p t denotes the average density of infected individuals (see 
Methods section for the numerical simulation method). Figure 2 
demonstrates that the prevalence declines with an increase in the 
low-risk patch ratio p in both homogeneous and heterogeneous 
metapopulations with the same mean degree. The homogeneous 
networks are given by the Erdos-Renyi random graphs 38 and the 
heterogeneous networks are by the Barabasi-Albert scale-free 
networks that are generated using a preferential attachment rule 39 . 
In the homogeneous networks, the decline rate of the prevalence of 
infection is higher than that in the heterogeneous ones. We can 
clearly observe an existence of the critical low- risk patch ratio p c , 
above which the prevalence is zero in the homogeneous networks. 
In contrast, the prevalence does not appear to fall to zero even for a 
sufficiently large p in the heterogeneous networks. The result 
suggests that the random intervention is not effective for the 
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Figure 1 | Metapopulation model with interventions for epidemic 
control in local patches. Schematic illustration of a metapopulation model 
consisting of high-risk patches without interventions (red open circles) 
and low-risk patches with interventions (blue open circles). The 
individuals in the patches are either susceptible (blue filled circles) or 
infected (red filled circles). In each patch, the disease transmission and 
recovery of individuals occur according to the SIS compartment model. 
The local basic reproduction number in the high-risk patches is given by 



i?< H) = fj K /n > l and that in the low-risk patches by B£> = < 1. The 
parameter p represents the ratio of the low-risk patches to the whole 
patches, corresponding to the total size of public health interventions. 

containment of epidemics if the patch connectivity is heterogeneous 
as in real transportation networks. 

For the random intervention, we analytically derive the critical 
low-risk patch ratio p c using heterogeneous mean-field approxi- 
mation 22,4042 . For the mean-field approximation to work effectively, 
we assume uncorrelated networks with a sufficiently large number N 
of the patches. When the mobility rates are positive, i.e., D s > 0 and 
Dj > 0, we obtain the critical low-risk patch ratio p c as follows (see 
Methods section and Supplementary Information for details): 
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Figure 2 | Transition between endemic and disease-free states for 
random intervention. The prevalence of infected individuals, Pi/p, is 
plotted against the low-risk patch ratio p. The circles and squares indicate 
the results for homogeneous and heterogeneous networks with N = 1000 
and (k) ~ 10, respectively. For each network type, the results for ten 
different network realizations are superimposed. The parameter values are 
set at p = 1, D s = 0.1, D, = 5, ft = 1, j8 H = 1.2, and p\ = 0.2. 



where P(k) represents the degree distribution and (k) = ■_ 1 ki/N 
represents the mean degree. A necessary condition for Eq. ( 1 ) is given 
by the inequality 

(k) n + D l 



P< 



(3) 



where fc max represents the maximum degree in the network. If 
inequality (3) does not hold, there exists a high-risk patch with 
degree k for which the denominator of the right-hand side of equa- 
tion (2) vanishes. In such a case, F(Pu) diverges, and thereby, p c goes 
to 1. The condition given by inequality (3) is equivalent to an exclu- 
sion of the possibility that the disease-free state is unstable for anyp 
in the range of 0 £ p < 1. See Supplementary Information for details 
and the results for the case with limited transmission. 

Figure 3 shows the phase diagrams in the (D^p) -parameter plane. 
The colors indicate the prevalence of infection in the steady state, 
separating the endemic (red) and disease-free (blue) states. The 
boundary between the two states corresponds to the critical low-risk 
patch ratio p c . The theoretical values ofp c obtained from equation (1) 
are in good agreement with the boundary in both the homogeneous 
and heterogeneous networks. 

Strategic interventions. Next, let us examine strategic interventions. 
The critical low-risk patch ratio p c depends on the procedure of 
selecting the low-risk patches with an increase in p. We consider 
three different types of intervention strategies that can reduce the 
value of p c as compared to the random intervention. An 
neighborhood intervention is inspired by an immunization strategy, 
called acquaintance immunization, on complex human contact 
networks 35 . In this strategy, all patches are assumed as high-risk 
patches initially at p = 0. Then, a patch is randomly chosen from 
the whole patches and one of its neighboring high-risk patches is 
changed to a low-risk patch at random. This procedure is repeated 
until the number of the low-risk patches reaches pN. A highly 
connected patch tends to be chosen as a low-risk patch by this 
procedure. The advantage of this intervention is that it only 
requires information about local patch connectivity. A targeted 
intervention is conducted based on centrality, which indicates the 
relative importance of the patch node from the viewpoint of network 
connectivity 43,44 . We only show the results for the usage of degree 
centrality, because the same results were obtained for the usage of 
closeness, betweenness, and eigenvector centralities. In this strategy, 
a high-risk patch with a higher centrality is more likely to be chosen 
as a low-risk patch. The targeted intervention requires information 
about global patch connectivity to compute the centrality values of 
the patches 33 . An adaptive intervention is introduced to test an ad- 
hoc intuitive manner in which interventions are targeted at unsafe 
patches with high morbidity rates. In this strategy, a high-risk patch 
with a high morbidity rate is preferentially chosen as a low- risk patch. 
Because morbidity rates in local patches vary with time, the choice of 
low-risk patches is adaptively updated at a fixed time interval T. 

Effectiveness of strategic interventions. We numerically examine 
the effectiveness of the neighborhood, targeted, and adaptive 
interventions in comparison to the baseline case with the random 
intervention. The best strategy, which enables us to contain 
epidemics with the smallest number of low-risk patches, can be 
identified in our simulations. The simulations were performed for 
different networks: homogeneous networks 38 , heterogeneous 
networks 39 , and an US airline network 19 . We introduce a parameter 
Aft to simultaneously change the transmission rates in both the high- 
risk and low-risk patches; they are given by p? H = 1 + AfS and fS L = 
Aft, respectively. As A/? increases, the transmission rate increases in 
all patches, and thereby make it difficult to contain an epidemic. 

Figure 4 shows the phase diagrams separating the endemic and 
disease-free states in the parameter plane of Aji andp for the homo- 



SCIENTIFIC REPORTS | 4:5522 | DOI: 1 0. 1 038/srep05522 



3 





10 u 
10" 
10" : 
10" : 
-j 10" 
10 _i 
10-' 
10"' 



Figure 3 | Validation of the analytical expression of the critical low-risk patch ratio for random intervention. Phase diagrams on the (D h p) parameter 
plane, where Dj is the mobility rate of the infected individuals and p is the low- risk patch ratio. The color bar indicates the numerically obtained values of 
the prevalence p x l p, averaged over ten different network realizations. The red and blue regions correspond to the endemic and disease-free states, 
respectively. The boundary between the two states corresponds to the critical low-risk patch ratio p c . The crosses indicate the theoretical values of p c 
obtained from equation (1) for the same ten network realizations. The parameter values are set at p = 1, D s = 0.1, /( = 1, /} H = 1.2, and fi L = 0.2. 
(a) Homogeneous networks with N = 1000 and (k) ~ 10. (b) Heterogeneous networks with N = 1000 and (k) ~ 30. 



geneous networks. It is reasonable that more interventions are neces- 
sary to contain the epidemic for a larger value of A/i in all the cases. 
Figure 4a shows the result for the random intervention, which is the 
baseline for comparing all intervention strategies. The result for the 
neighborhood intervention is almost similar to that for the random 
intervention, as shown in Fig. 4b. The targeted intervention based on 
the degree centrality is slightly better than the random intervention, 
as shown in Fig. 4c. Figures 4d, 4e, and 4f present results for the 
adaptive interventions with different update intervals T = 10, 50, and 
100, respectively. Overall, strategic interventions are not very bene- 
ficial for homogeneous networks with a narrow degree distribution. 

Figure 5 shows the phase diagrams for heterogeneous networks 
with broad degree distributions. As shown in Fig. 5a, the random 
intervention fails to contain the epidemic unless the interventions are 
given to nearly all patches. The neighborhood intervention is much 
better than the random intervention as shown in Fig. 5b. This is 
because high-degree patches containing a large population of indi- 
viduals tend to be chosen as low-risk patches in this strategy. 
Figure 5c shows that the targeted intervention based on the degree 
centrality is highly effective. The targeted interventions based on 
closeness, betweenness, and eigenvector centralities produced almost 
the same results, because all centrality measures are positively corre- 
lated to each other for the heterogeneous networks. Figures 5d, 5e, 
and 5f show the results for the adpaptive interventions with different 
update intervals T = 10, 50, and 100, respectively. These interven- 
tions are better than the neighborhood intervention, but not than the 
targeted intervention. In particular, the frequent update of low-risk 
patches in the case of T = 10 appears to be the best. The time 
evolutions of the average density of infected individuals exhibit oscil- 
lations due to the switching of dynamics, resulting from the update of 
low-risk patches (see Supplementary Fig. S3). The interplay between 
the continuous dynamics of epidemic spreading and the discrete 
switching of the low-risk patches causes a non-monotonic increase 
in the critical low-risk patch ratio with an increase in A/? as found in 
Figs. 5e and 5f. We can conclude that all the above mentioned stra- 
tegic interventions are highly effective compared to the random 
intervention, and the best strategy is the targeted intervention in 
the heterogeneous networks. 

To confirm the efficacy of strategic interventions in real hetero- 
geneous networks, we performed similar numerical experiments 
using the airline transportation networks in the United States 19 . 
The total number of airports is 500 and the total number of flight 
routes is 2980. The distribution of the number of flight routes shows a 



power-law decay, implying that the network has a scale-free prop- 
erty. Figure 6 shows the results for the random and strategic inter- 
ventions. The strategic interventions are considerably better than the 
random intervention, as observed in the previous results for the 
model networks with the scale-free property. However, in the US 
airline network, the adaptive interventions are less effective than 
the neighborhood intervention. 

Identification of target areas of intervention. Finally, we 
demonstrate that our framework is useful to identify the target 
areas for public health intervention using real human mobility 
patterns on the domestic airline network in Japan. From the time 
tables of domestic flights in All Nippon Airways (http://www.ana.co. 
jp) and Japan Airlines (http://www.jal.co.jp), we obtained the data 
about the total connectivity of the airline networks. Figure 7a 
illustrates the map of the airports in Japan. Further, we estimated 
the weights of the connections from the number of passengers 
traveling through each route in one day. The number of airports is 
N = 74 and the total number of routes is 190. The degree distribution 
is shown in Fig. 7b, approximately following a power law. The mean, 
maximum, and minimum degrees are given by (k) ~ 5.1, fc max = 50 
(corresponding to the Tokyo-Haneda airport), and k min = 1, 
respectively. Both the connection and weight matrices are 
symmetric. Because the weights are different between the flight 
routes as shown in Fig. 7b, we employ the metapopulation model 
with heterogeneous diffusion (see Methods section). The patches are 
assumed to correspond to the local cities around the airports. The 
initial density of individuals in each patch is set to be proportional to 
the total number of passengers at the corresponding airport. 

Figure 7d shows the prevalence of infection with the variation of 
the low-risk patch ratio p under the random and strategic interven- 
tions. The targeted intervention based on the degree centrality is the 
most successful strategy for the containment of epidemics by con- 
trolling only 5% of the patches, i.e., Tokyo (Haneda), Okinawa 
(Naha), Osaka (Itami), and Sapporo (New-Chitose). The second 
option would be the neighborhood intervention. The value of p c 
ranged between 8% and 16% for five simulations of this strategy. 
The adaptive intervention is not considerably different from the 
random intervention. Thus, our framework is useful to specify the 
best intervention strategy and the target areas where the control 
measures should be preferentially concentrated. If available 
resources and costs for epidemic control are limited, the best option 
can be found by comparing the prevalence for different interventions 
at a fixed value p. 
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Figure 4 | Effects of epidemic interventions in metapopulation models with homogeneous patch connectivity. Each panel shows the phase diagram for 
homogeneous networks in the (A/?, p) parameter plane, where A/J controls the transmission rates given by /? H = 1 + A/? and /? L = A/?. The color bar 
indicates the numerically obtained values of the prevalence Pi/p, averaged over ten different network realizations with N = 1000 and (fc) ~ 10. The 
red and blue regions correspond to the endemic and disease-free states, respectively. The parameter values are set at p = 1, D s = 0.1, and Dj = 5. 
(a) Random intervention; (b) Neighborhood intervention; (c) Targeted intervention; (d) Adaptive intervention with T = 10; (e) Adaptive intervention 
with T = 50; and (f) Adaptive intervention with T = 100. 
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Figure 5 | Effects of epidemic interventions in metapopulation models with heterogeneous patch connectivity. The same as Fig. 4, but for 
heterogeneous networks with N = 1000 and (k) ~ 30. The parameter values are set at p = 0.4, fi = 1, D$ = 0.1, and Dj = 5. 
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Figure 7 | Epidemic interventions in metapopulation models with the Japanese airline network. The data of the Japanese airport network were obtained 
from the timetables of domestic flights of All Nippon Airway (http://www.ana.co.jp) and Japan Airline (http://www.jal.co.jp). The Japanese airport 
network consists of JV= 74 airports and 190 flight routes. The mean, maximum, and minimum degrees are given by (k) — 5.1, fcn, M = 50 (corresponding 
to the Tokyo-Haneda airport), and fc m ; n = 1, respectively, (a) The map of the airports in Japan, which was created with Adobe Illustrator; (b) The degree 
distribution of the Japanese airport network; (c) The distribution of the number of passengers per flight route; and (d) The prevalence pj p is plotted 
against the low-risk patch ratio p for the random (squares), neighborhood (triangles), targeted (diamonds), and adaptive (circles) interventions. For the 
random and neighborhood interventions, the best result among five simulations was adopted. For adaptive intervention, the time interval for updating 
the low-risk patches is set at T = 10. The other parameter values are set at p = 0.5, ft = 1, D s = I\ = 10, /J H = 2, and /? L = 0.1. 
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Discussion 

We have presented a mathematical framework to evaluate the effec- 
tiveness of the random and targeted epidemic interventions in the 
SIS metapopulation models. First, we have obtained a theoretical 
result on the critical low-risk patch ratio for the random intervention 
and confirmed its validity by performing numerical simulations. We 
have shown that the heterogeneity of patch connectivity makes it 
difficult to contain epidemics. Then, we have compared the efficacy 
of intervention strategies using the result for the random interven- 
tion as a baseline. We have found that strategic interventions are 
considerably more effective than the random intervention in hetero- 
geneously connected networks. Moreover, the targeted intervention 
based on the degree centrality has shown to be most effective com- 
pared to the neighborhood and adaptive interventions. This result 
suggests that epidemic control measures should be focused on the 
highly connected areas. In contrast, the adaptive intervention can be 
worse than the neighborhood intervention. Our investigation has 
highlighted the importance of strategic interventions for epidemic 
control in mitigating and containing global epidemics. The optimum 
strategy can be selected with the support of the mathematical frame- 
work presented in this work. 

Our framework is applicable to more extended situations. The SIS 
model for disease transmission can be changed to other compart- 
ment models such as susceptible-infectious-recovered (SIR) and sus- 
ceptible-exposed-infectious-recovered (SEIR) models. For these 
models, an appropriate measure for the epidemic size should be 
the cumulative number of infected individuals during the period of 
an outbreak, instead of the prevalence that is used for the SIS model. 
We have considered the two transmission rates in the local patches 
for simplicity. A full exploration for the impact of the combination of 
the transmission rates on epidemic spreading is an issue to be studied 
in the future. Further, a more heterogeneous distribution of the 
transmission rates in the patches is a realistic factor to be incorpo- 
rated into the model. We have assumed a simple human mobility 
pattern (diffusion process) as a first step to test our framework. The 
human mobility pattern can be replaced by more realistic mobility 
patterns such as recurrent 15 , natural 27 , and safety-information-dri- 
ven 29 types. Furthermore, it is possible to consider other intervention 
strategies which have not been tested in this work; the preferential 
order of intervention can be determined based on other information 
characterizing the patches, such as the local population density. 

Methods 

Continuous-time metapopulation model. A metapopulation model consisting of N 
patches is considered. In each patch, the infection and recovery of individuals occur 
according to the susceptible-infectious-susceptible (SIS) model. In the continuous- 
time formulation for the progress of diseases on metapopulations, the equations for 
the dynamics in patch i (1 ^ i ^ N) are described as follows 24,25 : 

dp w 

-tt- = -PiPs,tPi,i + ViPu- D slh,i + D s^2 a P d jiPs,j> ( 4 ) 



=PiPs,iPi,i-PiPi,i-DiP 1J + D 1 '^2a ji d ji p 1 j, (5) 

where p S) , and p lti represent the numbers of susceptible and infected individuals in 
patch i, respectively. The epidemiological parameters, and p,-, denote the 
transmission and recovery rates in patch /', respectively. The term p'/Ps.iPi.i represents 
the nonlimited transmission, whereas the transmission term is given by PiPs,iPi,i/(ps,i 
+ PiJ) for the limited transmission 19,24,25 . The recovery rate is assumed to be the same 
in all patches, i.e., p,- — p. The transmission rate is set at = /? L < p for low-risk 
patches and at /?,- = fi H > p for high-risk patches. The ratio of low-risk patches to the 
whole patches is denoted by p. The parameters D s and D\ denote the mobility rates of 
susceptible and infected individuals, respectively. The bidirectional connectivity of 
the patches is characterized by the N X N adjacency matrix A = (fly): a,j = dp — 1 if 
patch i and patch j are connected; otherwise = a« = 0. The degree of patch /' is given 

EN 
,_ l a,j. When the diffusion is heterogeneous, connection weights are 

characterized by the N X N weight matrix W = (w«). The weight W//=Ww) is 
proportional to the number of passengers traveling between patch i and patch j. The 



diffusion rate is given by (i) dy — for the homogeneous diffusion and (ii) 
djj — Wjj j w $) ^ or t ^ le netero g eneous diffusion 23 . 

Numerical simulation method. We simulated equations (4)-(5) with the fourth- 
order Runge-Kutta method with time step 0.01 for f e [0, f max ], where f max = 2000. 

The total density of the individuals, p — . _ ; (p s f (f) + p l f (r)) ^N, remains 
constant over time t. The epidemic size was measured by the prevalence of infection 
Pi/p, where pj — Pi i( tmax ) / N is the average density of infected individuals for 

the random, neighborhood, and targeted interventions. For adaptive intervention 
producing oscillatory behavior of p I (f) with period 2T, the prevalence was measured 

byp I =(i/2T)£;:;-+ 2T £~ i p 1 , 1 ( t) /N. 

Mean-field approximation. As the low-risk patch ratio p increases from a sufficiently 
small value, the initially unstable disease-free equilibrium is stabilized at a critical 
value p — p c . Using the heterogeneous mean-field approximation 22 40 " 42 , we 
analytically derive the critical low-risk patch ratio p c for the metapopulation model 
(4)-(5) with uncorrelated networks and the homogeneous diffusion, i.e., (L — 1/fc*. 
We consider the case with positive mobility rates D s > 0 and Dj > 0 (see 
Supplementary Information for the other cases). We assume that inequality (3) is 
satisfied to exclude the case in which the disease-free equilibrium is unstable for anyp 
in the range of 0 < p < 1 . Given the average density of infected individuals in the 

steady state, p Y — ._ i p Y j^N, we can rewrite the model (4)-(5) as follows: 

~^ L = -A^s,iPu + m, I -- D sPs,,+^si|r(^-Pi), (6) 



-rf = PiPsjPu - m,i - APi,i + A 77T Pi- ( 7 ) 



The endemic equilibrium of the above equations is given by the positive solution of 
the following equation: 

where D 1S = Di/D s and y(Pi) = p + (D r s — l)p : . See Supplementary Information for 
the explicit expression of p\ f as a function of p { . The self-consistency condition for p t 
yields the following equation: 

Pi = 4(l>M(A) + I>y(ft)). W 

where H and L denote the sets of indices for the high-risk and low-risk patches, 
respectively. Denoting the right-hand side of equation (9) by G(pi), the condition for 
the critical point is given by G'(0) — 1. We can evaluate G'{0) as follows: 

G(0)~(l-p)F(ll H )+pF(li L ), (10) 

where F{f3) is denned in equation (2). By solving G'(0) — 1 with respect to p, we obtain 
the formula for the critical low-risk patch ratio given by equation (1). When 
inequality (3) does not hold, F(/? H ) diverges, and hence, the critical value p c goes to 1. 
See Supplementary Information for the analysis of the case with limited transmission. 
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